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dominates. We compute the complete set of 0{as) SQCD corrections to both production channels 
and study their dependence on the parameters of the Minimal Supersymmetric Standard Model. 
In particular, we discuss the prospects for observing strong, loop-induced SUSY effects in top-pair 
production at the Tevatron Run II and the LHC. 
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I. INTRODUCTION 

Since the discovery of the top quark at the Fermilab Tevatron pp coHider in 1995 
the top-physics program has shifted to precisely studying its properties. The high top- 
quark yields at the Tevatron Run II and the CERN LHC pp collider open a new, rich 
field of top-quark phenomenology, which may enable a precision physics program with top- 
quark observables, in the spirit of the successful LEP/SLC studies of precision Z boson 
observables j^. 

The Standard Model (SM) of electroweak and strong interactions has seen impressive 
experimental confirmation. However, the Higgs boson, which is predicted by the SM as 
a direct consequence of mass generation via spontaneous electroweak symmetry breaking 
(ESB), has so far eluded direct observation. Moreover, the many shortcomings of the SM 
suggest that it may indeed be a low-energy limit of a more fundamental theory. One of the 
most promising candidates for a theory beyond the SM is Supersymmetry (SUSY) (for a 

nnnpn 

review see, e.g., Refs. [J, 15|, lo, lZ|, 18||). Supersymmetry solves the finetuning problem, allows 
for gauge coupling unification, provides a dark matter candidate, predicts a light Higgs 
boson and agrees with precision electroweak measurements. However, no direct experimental 
evidence of SUSY has been found yet. 

The large mass of the top quark suggests that it plays a special role in the mechanism of 
ESB, and that new physics connected to ESB may be found first through precision studies 
of top-quark observables. Deviations of experimental measurements from the SM predic- 
tions, including electroweak and QCD corrections, could indicate new non-standard top 
production or decay mechanisms. If supersymmetric particles are detected at the LHC, 
the comparison of precisely measured top-quark observables with their predictions including 
SUSY loop effects may yield additional information about the underlying model, that may 
not be accessible otherwise. Measuring precisely the properties of the top quark therefore 
is an important goal at the Tevatron Run II and the LHC. To fully exploit the potential 
of these colliders for precision top-quark physics, it is crucial that predictions for top-quark 
observables include higher-order corrections within the SM and beyond. In this paper we 
study the impact of SUSY QCD (SQCD) one-loop corrections within the Minimal Super- 
symmetric SM (MSSM) on strong top-pair production at both the Tevatron Run II and the 
LHC. 
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Presently, the total top-pair production cross section has been measured at the Tevatron ^ 
with a relative uncertainty of Aatj/crtf = 12 % ( with L = 760 pb^^) and is in good 



agreement with the theoretical QCD prediction [10|,llll|- It is anticipated that at the Tevatron 
Run II the top-pair production cross section will be measured with a relative uncertainty of 
Aatt/cytt ~ 10% (with C = 2 fb~^). At the LHC, the goal is to measure aa ultimately with 
a relative uncertainty of < 5%. The current experimental uncertainty still leaves room for 
SUSY loop effects. Through the virtual presence of supersymmetric particles in quantum 
corrections, the measurement of the cross section, kinematic distributions and the extraction 
of masses and couplings may be affected. Among the top-quark observables under study, 
parity-violating asymmetries in the production of left and right-handed top quarks are of 
special interest. They have the potential to provide a clean signal of non-SM physics: QCD 
preserves parity and the SM induced asymmetries are too small to be observable, at least 
at the Tevatron pp collider l^- The produced top quarks decay almost entirely into 
a bottom quark and a W boson before they can hadronize jl^ or flip their spins. The 
spin correlation of the top-pair system will therefore be preserved and can be measured 
by studying angular distributions of the decay products [l5|. In order to determine the 
sensitivity of top-quark observables to SUSY loop effects, it is necessary to calculate the 
radiative corrections to the top production and decay processes, both within the SM and its 
supersymmetric extension, and to implement these calculations in a Monte Carlo program. 
The latter will allow an efficient determination of those observables that are most sensitive 
to supersymmetric particles, after taking into account the detector response. 

In this paper, we present the calculation of the complete SQCD 0{as) corrections to 
the main production channels of strong top-pair production at the Tevatron and LHC, 
qq — > tt (qq annihilation) and gg ti (gluon fusion), and present a detailed study of their 
numerical impact on top-quark observables in both unpolarized and polarized ti production 
(the initial-state quarks and gluons are unpolarized). We consider the MSSM with CP- 
conserving couplings. The inclusion of the top-quark decays and implementation in a Monte 
Carlo program for off-shell ti production is work in progress {l^, but beyond the scope of 
this paper. 



^ See www-cdf.fnal.gov and www-dO.fnal.gov for most recent results from the CDF and D0 collaborations, 
respectively. 
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The SM predictions to both polarized and unpolarized tt production at hadron coUiders 
are under good theoretical control: The Born-level amplitudes to qq annihilation and gluon 
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22|. At next-to-leading order (NLO), 



fusion were first considered in Refs. Q, 

the SM QCD corrections for the total cross section and s pin -independent kinematic distri- 

28|, electroweak corrections have 



been studied in 


Refs. 
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effects in Refs. 
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35l | and soft gluon resummation and threshold 
To measure spin correlations and asymmetries at the 



Tevatron and LHC, higher-order corrections to polarized tt production need to be known 
as well: The NLO QCD and electroweak corrections for the ^larized production have been 

and Refs. 



calculated in Refs. 



40 



41 



131, 



30 



32 



4j|, respectively. 



A number of studies of higher-order corrections within the MSSM are also available: 
The SUSY electroweak one-loop corrections to unpolarized on-shell ti production at hadron 



colliders have been calculated in Refs. 



45 



46 



47| {qq annihilation) and Ref. [47| (gluon 



fusion), and the same set of corrections for polarized top-pair production has been considered 
in Ref. |12|. A study of leading logarithmic SUSY electroweak corrections can be found 
in Ref j48| and of Yukawa corrections in a two-Higgs-doublet model is given in Ref. j^^. 
The NLO SQCD corrections to the gg ti subprocess without polarization have been 
presented in Refs. l5ll | and in Ref. [511 for the scattering of polarized protons. The 



unpolarized qq annihilation subprocess has been investigated at NLO SQCD by a number of 
groups: In Ref. H the corrections with the assumption of negligible box contributions have 



53| the gluon self-energy and crossed box diagram 



been calculated, whereas in Refs. 

contributions have been neglected. The full NLO SQCD corrections to qq tt have been 



examined in Refs. 



However, the results of these difi'erent calculations do not 



agree. The origin of this disagreement has been found (see footnote in Ref. and is due 

to a difference in the signs of the direct and crossed box diagrams (see Fig. [4]), leading to 
large numerical differences that are especially pronounced in case of comparable small top 
squark and gluino masses. However, a reliable, complete and detailed study of the impact 
of NLO SQCD corrections on ti production based on the correct result is not available in 
the literature. 

This paper represents the first detailed study of the numerical impact of NLO SQCD 
corrections on both polarized and unpolarized strong ti production at hadron colliders, that 
considers both main production channels, qq annihilation and gluon fusion. We present nu- 



4 



merical results for the total production cross section, the invariant ti mass and top transverse 
momentum distributions, and polarization asymmetries, at both the Tevatron Run II and 
the LHC, respecting current experimental bounds on the masses of the supersymmetric par- 
ticles. We provide analytic expressions for the complete NLO SQCD corrections to polarized 
ti production and compare our results with the literature where available. Our results have 
been thoroughly cross-checked by performing at least two independent calculations, based 



on analytic methods and the FeynArts 



57|, 



58| and FormCalc 



59|, 



GOfl packages. 



The paper is structured as follows: In Section [TTl we present the analytic results for 
the partonic qq, gg ti scattering processes, specify our renormalization procedure and 
our convention of the squark mixing parameters and choice of MSSM input parameters. 
Section UTTl contains a detailed study of the numerical impact of the NLO SQCD corrections 
on polarized and unpolarized top-pair production at the Tevatron RUN II and LHC. We 
first discuss their impact in detail at parton level and then present their most pronounced 
effects to the observable, hadronic cross sections, followed by a discussion of their prospects 
to be observable at these colliders. We also include a comparison with available results in 
the literature. Finally, a summary of our results and conclusion can be found in Section HVl 
Details of the calculation and explicit analytic expressions are provided in the appendices. 
We also made available a Fortran code ^ that calculates the cross sections presented in this 
paper to polarized and unpolarized top-pair production at hadron colliders at leading-order 
(LO) QCD and NLO SQCD. 

II. STRONG TOP-PAIR PRODUCTION AT NLO SQCD 

At the Tevatron the main production mechanism for the strong production of top-quark 
pairs is the annihilation of a quark-antiquark pair 

liP-i) + q{P3) 't{p2, Xt) + i{pi, Xt) , 

whereas at the LHC the top-quark pairs are mainly produced via the fusion of two gluons 

giPi) + gips) ^t{p2, Xt) + i{pi, Xt) , 



The Fortran code is provided at http://ubpheno. physics. buffalo. edu/'-^^dow/ppttsqcd. 
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where pj, i = 1, ■ ■ ■ , 4 are the four-momenta of the in- and outcoming particles and Af(Aj) = 
±1/2 denotes the top(antitop) hehcity state. The corresponding Feynman diagrams at LO 
QCD are shown in Figs. [H [2l The partonic differential cross sections to the qq annihilation 
and gluon fusion processes for polarized top-quark pairs at NLO SQCD can be written as 



m,99 



"'^gg,99 



(t, s, At, Af) + ddaqq^gg{i, s, Xt, Af) 
"^IMfl^ + 27^e ^(^A^J^S'' X Mf'''*' 



(1) 



where s = {pi + P2) = (ps + Pi)\t = {ps - piY = (p4 - P2) = mf - s{l - /3t cos6)/2 
are Mandelstam variables with 6 denoting the scattering angle in the parton center-of-mass 
system (CMS) and fit = a/1 — 4m^/s is the top-quark velocity. The matrix elements squared 
are averaged over initial-state spin and color degrees of freedom and summed over final-state 
color degrees of freedom. The phase space of the 2 — > 2 scattering process, ^$2^2, as usual 
reads 



2^2 



'{P3+P4 



Pi - P2) 



Pt 



2n 



d cos 6 



(2) 



2p? 2p0 

where the phase space has been evaluated in the parton CMS with 0* denoting the azimuthal 
angle. The matrix elements, and SAdg^^^g^ , describe respectively the LO QCD and 

one-loop SQCD contributions to qq annihilation and gluon fusion. The explicit expressions 



for M''^'^^ anc 
found in Ref. 



the spin and color averaged transition amplitude squared, XII-^b'^^P? can be 



29|| (Sec. 2.1 and 2.2) for unpolarized and in (Appendix A) for polarized 



top-pair production. 

The observable hadronic differential cross sections are obtained by convoluting the par- 
tonic cross sections of Eq. ([T]) with parton distribution functions (PDFs) 



da 



LO,NLO 



1 + 



ij=qq,99 



ij Jo 



dxidx2 



fi{xi,^F)fj{x2,fiF)dcr^^ 



LO,NLO 
ij 



(asifiR), s, i, At, At) + i ^ J 



(3) 



with S = s/{xiX2). In the numerical evaluation, we use the CTEQ6L1 (LO) set of parton 



distribution functions 



6l| with the QCD factorization [jj^f) and renormalization scale (jj^r) 



chosen to he fiF = fiR = mt. Accordingly, the cross section is evaluated using the one-loop 
evolution of the strong coupling constant with nj = 5 light flavors, ^qcd — 0.165 GeV and 



af°(Mz) = 0.130, which yields af^(mt) = 0.1176. 
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In the remaining parts of this section we will present the explicit expressions for the par- 
tonic NLO SQCD differential cross sections, da^^^^ and da^J"'^, to polarized top-quark pair 
production. The corresponding cross sections for unpolarized ti production can be obtained 
by summing over the top and antitop-quark helicity states, J2xt At=±i/2 d^q^gg'^^i^t, Af). 




Figure 1: Feynman diagram to the qq tt subprocess at LO QCD. 




(a) (b) (c) 

Figure 2: Feynman diagrams to the s (diagram (a)), t (diagram (b)) and u (diagram(c)) channels 
of the gg — > tt subprocess at LO QCD. 



A. NLO SQCD corrections to qq annihilation and gluon fusion 

The SQCD 0{as) corrections modify the tree-level gtt,ggg and gqq vertices and the 
gluon propagator through the virtual presence of gluinos (^), squarks {qL,R), stops (tL,i?) 
and sbottoms {Pl,r)^ i-e. the superpartners of the gluon, the left and right-handed light 
quarks, top and bottom quarks, respectively, as shown in Figs. [SMl For completeness, we 



also provide the corresponding 



We closely follow Refs. 



12 



bynman rules in Appendix [Al 



29| and write the NLO SQCD matrix elements, dM^^^^ 



'■qq,gg ' 



of Eq. dl]) in terms of form factors that describe the SQCD one-loop corrections. In case 
of gluon fusion, these form factors multiply so-called standard matrix elements (SMEs) 
that contain the information about the Dirac matrix structure. After the interference with 
the Born matrix elements, the SMEs are written in terms of scalar products involving the 
external four-momenta and the top/antitop spin four- vectors, s^^. The latter are defined 
after choosing the axes along which the t and i spins are decomposed, as described in 






(a) (b) (c) 

Figure 3: Generic self-energy (diagram (a)) and vertex (diagrams (b) and (c)) corrections to qq — > ti 
at NLO SQCD. The corrections to the gluon propagator and the gqq vertex are explicitly shown in 
Figs. [71 and\27\ respectively. 

1. . ^.t 




I 9 , 

^" . : ^ : 

(a) (b) 

Figure 4: (a) Direct box diagram and (b) crossed box diagram contributing to qq tt at NLO 
SQCD. Graphs containing squarks/stops are summed over the squark/stop mass eigenstates i,j=L,R 
(no mixing), i,j=l,2 (with mixing). 



Appendix A of Ref. [12||. As studied in Ref. [43||, for instance, the freedom in the choice of 
the spin axes can be used to increase spin correlations at hadron colliders. Here we choose 
the helicity basis where the spin is quantized along the particle's direction of motion. Using 
the helicity basis, the NLO SQCD contribution to the polarized partonic cross sections of 
Eq. ([T]) read as follows (with z = cos 9): 
qq annihilation: 

27^e '^{6M^g9C'' X Mf) = 

E I I' (^Fv{s,m, = 0) + Fv{s,mt) - fl{s)^ + 

— ^7^e - z^)il + A\tXt)FM{s,mt) - 2(At - Af) [2zGA{s,mg = 0) + 
y 

where the gti and gqq vertex and the gluon self-energy corrections shown in Fig. [3] are 
parametrized respectively in terms of UV finite (after renormalization) form factors Fy, 
Fm, Ga, and the subtracted gluon vacuum polarization n(s) = n(s) — 11(0). Fy, Fm, and 
Ga denote the vector, magnetic and axial vector parts of the vertex correction. The direct 
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(b) 




I crossed 
diagrams 






Figure 5: Generic vertex and self-energy corrections to the t{u) channel (diagrams (a)-(c)) and the 
s channel (lower row) of the gg ti subprocess at NLO SQCD. The corrections to the gti vertex, 
the quark and gluon propagators, and the ggg vertex are explicitly shown in Figs. [2% [M, fTj. and\2^ 
respectively. 

and crossed box contributions to qq annihilation shown in Fig. [4] are described by Bt and 
respectively. Explicit expressions for these form factors and the box contributions are 
provided in Appendix [B] and n(s) is given in Eq. ffT2l) . 
gluon fusion: 



27^e Y.i^M'g^'''' X M'i*) 



64 



2 7^e<{ (c^(j). 



i=l,...,7 
11,. ..,17 



..v,t .^ ( P^ii s) pf\i, s) > A,t 



<rr-'((,.5) , <Tf-'(i,j) 



i — m\ ' {i — m^Y 



V,u 

n — ml {u — mf^'^ 



V.u / ^ ^\ Tj,u / ^ ^\ 



+ 



u — mi {u — mi) 



2\2 



d=a,b,c 



d=a,b 



+ 



^[m1;1^_(12,2) + M,^1^_(12,3) 



(5) 
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9 

9 : 



(a) 



I 

(b) 



-|- crossed diagrams 




(c) 





(d) 



Figure 6: Box corrections to gg ti at NLO SQCD. The box diagrams d = a,b contribute to both the 
t and u(crossed) channels, whereas the box diagrams d = c and the box-triangle correction (diagram 
(d)) involving a quartic ggtt interaction only contribute to the t channel. Graphs containing stops 
are summed over the stop mass eigenstates j=L,R (no mixing), j=l,2 (with mixing) 



where i numerates the 14 SMEs, M 



iV,A),it,n) 



of Eqs. (Bl), (B2) in Ref. [29|, andj = 1,2, 3 the 



s, t, u channel of the Born matrix element to the gluon-fusion subprocess. For convenience, 
the color factors c'*'*'"(j) of Ref. 29|| and c^'"(j) are both provided in Appendix IC 31 The 



coefficients to the SMEs describe the parity conserving {Pi^'^^'^^'^\ p^d^''^\ Pu^^^^) and parity 



violating (cr. 



{v,s),(t,«) ^n,{t,u) 



parts of the SQCD one-loop corrections to the t and u channels 



of the gluon-fusion subprocess, which consist of 



• vertex corrections that modify the gtt vertex (p]^'^*'"-*, 0"^^'*^*'"''), shown in Fig. [5] (dia- 
grams (a) and (b)), 

• top-quark self-energy corrections that modify the top-quark propagator 
(pf'^*''^'*, erf ''■*'"''), shown in Fig. [5] (diagram (c)), 

• box contributions {p^d^'^\ o'^d*'^^), shown in Fig. [6] (diagrams (a), (b) and (c)), and 

• the box-triangle correction (p^^'^^**), shown in Fig. [6] (diagram (d)). 

The SQCD one-loop corrections to the s channel of the gg — > ti subprocess, shown in 
Fig. E] (lower row), are parametrized in terms of the gtt vertex form factors {Fv,Fm,Ga), 
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the coefficient describing the ggg vertex corrections (P2'^)i the subtracted gluon self- 



procedure. Their explicit expressions are provided in Appendix O The expressions for the 



Before we turn to the discussion of the numerical impact of these higher-order SQCD 
corrections on ti cross sections, we describe in detail in the following sections our choice of 
the renormalization scheme and the MSSM input parameters. 

B. Renormalization scheme 

The complete set of SQCD one-loop corrections to qq annihilation and gluon fusion is 
gauge invariant and IR finite. However, the self-energy and vertex corrections shown in 
Figs. [8f29l exhibit UV divergences that arise in form ofA = 2/e — 7^; + log(47r) terms when 
using dimensional regularization in = 4 — e dimensions. These singularities are removed 
by introducing a suitable set of counterterms that are fixed by a set of renormalization con- 
ditions. We employ multiplicative renormalization and perform the following replacements 
of the left and right-handed quark fields, \E'l,/? = |(1 =F 75) the top-quark mass, nit, the 
gluon field, GJJ, and the strong coupling constant, gs, in the QCD Lagrangian: 



With Zi = 1 + SZi, this yields the counterterms to the quark and gluon self-energies, qqg 
and ggg vertices as follows (with 6qt = for g 7^ t and a, b, c denoting respectively color 
degrees of freedom of quarks and gluons): 



energy (n(s)). All coefficients to the SMEs are UV finite after applying the renormalization 




(z, j), can be found in 





iSij^iSZy - SZa'J^) - 5qt{mq5Zv - Sniq)] 



(6) 





(7) 



k 





(8) 
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: -9s fabc [{ki - A;2)M3fi'/iiM2 



(9) 



with 5Zi = 5Zg + 35Z:i/2 and SZy^ = {5Zl±5Zji) /2. For the renormalization of the strong 
coupling constant and the g 



62 



uon field, we use the MS scheme, modified to decouple the 



heavy SUSY particles 
are subtracted at zero momentum: 

• oZ% = ■ 



63|, i.e. divergences associated with the squark and gluino loops 



-n(o) 



(10) 



Qf^2 \e,modified 

where TiT\e,modified dcuotcs the terms proportional to A^, = A — ln(m^//i|,) (with 
m = niq, nig) of the transverse part of the gluon self-energy, yielding 



2Ar. 



(11) 



Only the transverse part of the gluon self-energy contributes to the NLO SQCD cross 
section for ti production and the SQCD one-loop contribution to the renormalized 
gluon self-energy shown in Fig. [7] enters the partonic differential cross sections of 
Eqs. (HD, ([5]) as follows: 



as I \ ^ \ ^ 1 

47r I ^ ^ 6 

q — u,c,t j 
d,s,b 



(fc^ — 4m 



Bo{k , rrig, nig) + -p-5o(0, m^^., m^-J + - 



4m? 



2 (2nii + k^) „ 

3 + 2^ p -Bo{k\ m-g, m-g) - -j/Bo{0, m-g, m-g) \ - n(0) 



(12) 



where j = L, i?(l,2) (no(with) mixing) sums over the two 
and -Bo denotes the scalar two-point function (see, e.g., Ref. 



squark mass eigenstates 



SZi — —T\^^modified = ~n(0) , (13) 

where T\^^rnodified dcuotcs the term proportional to Am (with m = niq^nig) of the one- 
loop ggg vertex correction shown in Fig. [29] that multiplies the tree level ggg vertex. 
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Using Eqs. (fTOll . (fTSll . one can verify that SZi = SZ3 and, thus, SZg = —SZ3/2. It is 



interesting to note that in this modified version of the MS scheme the counter term for 



29l |. Moreover, there is no contribution from 



the qqg vertex coincides with the one of Ref. 
self-energy insertions into external gluon legs. 

To fix the renormalization constants for the external quarkSj_5Zy A(?Tiq), and the top-quark 



mass renormalization constant, 6mt, we closely follow Ref. [29| and use on- shell renormal- 
ization conditions. The renormalization constants (5ZyA(^g), ^m^ are determined by the 
vector, axialvector and scalar parts of the quark self-energy, evaluated at the on-shell quark 



mass, T.v,A,s{P^ 



for the top quark and T,v,a{p'^ 



0) for the incoming quarks in 



qq annihilation. The SQCD one-loop corrections modify these quark self-energy contribu- 
tions as shown in Fig. [8] as follows: 



Atx 



3 rrir, 



(14) 



yielding the following on-shell renormalization constants: 

d 



SZv(,m„ 



5ZA{mq 



Suit 



47r ^ 3 

j 

2 



(15) 



mt 



m„ 



^sip' 



_as_ \ - 2 

47r ^ 3 

i 



—X\Bi + A • — -Bq 
■' nit 



{nil,nig,nii 



(16) 
(17) 



where j = L,_R(1,2) (no(with) mixing) sums over the two squark mass eigenstates, and 
the coupling constant factors, Xf,Xf, are given in Appendix [Al The scalar and vector 



two-point functions, i?o,i, and their derivatives, B'q^{p^ 



ni 



found in Ref. 



dBo,i/dp^\p2 



64| . for instance. In the on-s/iel/ renormalization scheme, the derivative of the 



renormalized quark self-energy of Fig. [8] is zero when evaluated at = m^, and thus there 
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is no contribution from the self-energy insertions in the external quark lines. Finally, it is 
interesting to note that for our choice of renormalization scheme, the renormalized SQCD 
one-loop corrections only depend on through the strong coupling constant, as{fJ.R). 





9 b „a b a I \ b „a ,h 

9^9' 9 ^ g \ , 9 , 9 \/ 9 



Figure 7: Renormalized gluon self-energy , — lYj^i, — i^abif^fj.^v T^L-g^u^x), atNLOSQCD. Graphs 
containing squarks are summed over the squark mass eigenstates j=L,R (no mixing), j=l,2 (with 
mixing) and the quark flavors q={u,d,s,c,b,t}. 




Figure 8: Renormalized quark self-energy, iT, = — + mqT,s), at NLO SQCD. Graphs 

containing squarks are summed over the squark mass eigenstates j=L,R (no mixing), j=l,2 (with 
mixing). 



C. Scalar quark sector of the MSSM 

The superpartners of the left and right-handed top quarks are the left and right-handed 
scalar top squarks ('stops'), ii and i^. Assuming all supersymmetric breaking parameters 
are real, the part of the MSSM Lagrangian that contains the stop mass terms reads 

^= -iil,i*n)^ \ ~ ] ^ith M=\ (18) 

where the diagonal entries m?^ and m?^ of the stop mass matrix A4 as well as Xf are defined 
as 

ml = Ml+ m\ + M| cos 2/5 {l\ - Qtsl) , Xt = (At - fi/ tan (3) , 

ml = Ml+ ml + M| cos 2/5 Q^sl (19) 
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with /g = 1/2, Qt = 2/3 and = 1 — M^/M|. Furthermore, tan/? = t'2/t'i (also Vu/vd in 
the literature) denotes the ratio of the two Higgs vacuum expectation values, the trilinear 
Higgs-stop-stop coupling and fi is the Higgs-mixing parameter. M^^ and M^^ are the soft 
supersymmetry-breaking squark-mass parameters. As can be seen from Eq. (ITSl) . ii and 
iji are not necessarily mass eigenstates, since ^A is of non-diagonal form. Thus, ^l,^/? can 
mix, so that the physical mass eigenstates ti,t2 are model-dependent linear combinations 
of these states. The latter are obtained by diagonalizing the mass matrix performing the 
transformation 







[hi 





cos 6i sin 6f: 
— sin 9^ cos 9^ 



)l 






I in ) 



(20) 



The mass term of the supersymmetric Lagrangian will then transform to the diagonal form 

^ = -iii, I ""'^ \ I r: I • (21) 



m? 

Inserting Eq. (l20l) and its adjoint form into Eq. (l2Til and comparing with Eq. (fTHll the mass 
matrix M. can be written as 

\^ sin cos (m| - m| ) sin^ 9i ml + cos^ m| ^ 
Comparing the two different descriptions of the mass matrix we find 



2 2 /) 2 I • 2 /) 2 

It = cos mr + sm 9z rrir 



2 • 2 /I 2 I 2 /I 2 

rrir = sm t^fm,- + cos ^fm,- , 

nitXt = sin 9^ cos 9 f (ml — fn^) . (23) 

As input parameters we choose the physical measurable quantities m^^, mi^ and the stop 
mixing angle Studying the properties of Eq. (|23|) we find that the three equations, 
and therefore the cross sections of Eq. ([1]), are invariant under the following two sets of 
transformations: 

9i 9i + n-7i (nel) (24) 
m| ^ m| and 9^ ^ 9i + tx/2 . (25) 

Therefore, for a complete exploration of the size of the SQCD one-loop corrections, it is 
equivalent to either vary the masses m?^ and m?^ completely independently and vary the 
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mixing angle between ±'/r/4, or to choose one of the squark masses to be always the lighter 
one and vary 6i between ±7r/2. We choose the latter option and take m| always to be the 
lighter top squark and vary the stop mixing angle in the range — 7r/2 < 6i < 7r/2. 

Taking m^j, m^^ ^ind 9^ as input parameters and using Eq. (i23l ). we can calculate m?^, m?^ 
and Xt- The soft supersymmetry-breaking parameters Mi^ and Mi^ can be determined from 
Eq. ffTOll in dependence of tan /3. Since the SQCD NLO cross sections do not directly depend 
on tan [3 and the dependence of m?^ and m?^ on tan [3 in the range 1 < tan /3 < 50 is rather 
weak, we choose in this paper an arbitrary value for tan/5, i.e. tan/? = 10. 

The equations (fT8l)-(i23l) can be formulated in a similar way for the other squark flavors, 
one just has to replace the index "t" by the appropriate squark index, as well as Xt — * Xd = 
{Ad — fi tan P) for the down-type squarks. Since the off-diagonal elements of the squark mass 
matrices, m^Xq, are proportional to the quark masses m^, and Xq cannot be arbitrarily large 
{Xq < 65|), mixing can be neglected for the first two generations of squarks and may be only 
important in the sbottom sector for large values of tan /5. However, since the influence of 
the sbottom sector on the hadronic ti cross sections is very small, we also neglect mixing for 
the bottom squarks and in addition assume m^^ = mi^. Because gauge invariance requires 
M-t^ = , it follows that with 

ml =Ml+ml + Ml cos 2(3 {ll - Q,sl) , (26) 

/| = —1/2 and Qt, = —1/3, the bottom squark masses are already fixed by the input 
parameters of the top squark sector. We further assume universal squark masses for the 
first two generations of squarks: 

rriu^ = rriu^ = = = m^r. = ^cr = ^Sl = ^Sr ■= rriq . (27) 

To summarize, in this paper we will study the dependence of the NLO SQCD predictions 
on the following MSSM parameters: rrig, m^-^, m^-^, 9^ and rrig. 

III. NUMERICAL RESULTS 

In this section, we examine the numerical impact of the SQCD one-loop corrections on 
the cross sections for unpolarized and polarized top-pair production at the Tevatron Run II 
(v^ = 1.96 TeV) and LHC (v^ = 14 TeV). The details of the calculation are presented 
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in Section [TTl and the Feynman rules and explicit analytic expressions are provided in Ap- 
pendix [A] and [Bl O, respectively. 

We performed a number of checks of our calculation: The results for both the qq annihi- 
lation and gluon fusion matrix elements for polarized and unpolarized top-pair production 
have been derived by at least two independent analytic calculations based on the Feynman 
rules of Ref . 



66|. The numerical 



plementation described in Ref. 



evaluation of the tensor integrals was performed with an im- 



29| as well as by using the packages LoopTools 



6fl/FF 



67|. 



681 . The results 



Furthermore, one of the analytic calculations used the package TRACER 
of the analytic calculations have been compared numerically with the results obtained with 



the FeynArts 



57,3 



and FormCalc 



59 



60|| packages using the MSSM model file of Ref. |69| 



In all cases, we found perfect numerical agreement. 

In the computation of the ti observables, we use the following values for the SM input 
parameters [8|: 



mt = 175 GeV, rrib = 4.7 GeV, Mw = 80.425 GeV, Mz = 91.1876 GeV . 



(28) 



As described in Section [TTl (Eq. Q), we use CTEQ6L1 PDFs [GjJ and the strong coupling 
constant is evaluated at LO QCD with the factorization and renormalization scales chosen 
to be equal to the top-quark mass. Our choice of MSSM input parameters is described in 
Section III Ci As an additional constraint on the MSSM parameter space, we calculate the 
SUSY loop corrections to the p parameter, Ap{mi, rrii, 6^) (for a review see, e.g., Ref. (toI). 
and only allow for MSSM input parameters that yield Ap < 0.0035 [TjJ. If not stated 
otherwise, we choose the masses of the supersymmetric partners of the light quarks to be 



rrin 



2 TeV. In varying the remaining MSSM input parameters we observe the following 



experimentally motivated mass limits 



72l, 0, Q, l75fl: mg> 230 GeV and m^-^ > 100 GeV. 



In the following sections, we study the numerical impact of the SQCD one-loop corrections 
on the total ti production rate, the invariant ti mass and top transverse momentum distri- 
butions, and on polarization asymmetries. We discuss the case of unpolarized top quarks 
in Section IIII Al and consider polarized ti production in Section IIII B[ Before we present 
results for the observable hadronic cross sections at the Tevatron Run II and the LHC in 
Sections IlII A 31 ITlI B 31 we perform a detailed investigation of the impact of the SQCD one- 
loop corrections on the parton-level cross sections in Sections IlII A 11 [Til B l\ (qq annihilation) 
and Sections IlII A 21 [Tll B 21 (gluon fusion). Finally, in Section [Tll CI we compare our results 
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with existing calculations in the literature where available. 



A. Unpolarized top-pair production at NLO SQCD 

We first consider the produced top quarks to be unpolarized, i.e. the corresponding 
partonic and hadronic cross sections are obtained from the ones of Eqs. ((I|),((3l) by sum- 
ming over the top and antitop-quark helicity states, J2xt Xt=±i/2'^^q^gg^^'^i^^ -^t^ -^t) 
,At=±i/2 daLo,NLo{S, Xt, k), respectively. 

In order to reveal the numerical impact of the SQCD one-loop corrections on observ- 
ables to unpolarized top-pair production and to study the dependences on the MSSM input 
parameters, the gluino mass {rrig), the two top-squark masses {mi^^), and the stop mixing 
angle (6'j), we introduce relative corrections as follows: 

• At the parton level, we use 

A.,„,(^) = '^'^.o ' (29) 

where a^q^g^^^ {s) denote the total partonic production rates to unpolarized ti produc- 
tion via qq annihilation and gluon fusion at LO QCD and NLO SQCD, respectively. 

• At the hadron level, we compute the total hadronic cross section, the invariant ti mass 
{Mtt) and top transverse momentum [px) distributions at the Tevatron Run II and 
the LHC, and use 



tt) 



daio/dMu 



^KPt) = 1 r; , (30) 

daio/dpT 

where daLo,NLO describe the hadronic cross sections to unpolarized tt production at 
LO QCD and NLO SQCD, respectively. All hadron level results are obtained by 
including both processes qq ti and gg — > ti. 
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1. Effects of NLO SQCD corrections in qq annihilation for unpolarized top quarks 

At the Tevatron, top-quark pairs are mainly produced via qq annihilation and, thus, 
some of the characteristics of the corrections discussed here at the parton level will manifest 
themselves again in the hadron-level results for the Tevatron. In the following we study the 
impact of the NLO SQCD corrections on the uu ti subprocess, representatively for all qq 
initiated processes. 

Fig. [9](a) illustrates the impact of the (renormalized) gluon self-energy, shown in Fig. [7] 
(see also Eq. ffT2l)). in dependence of the partonic CMS energy, The upper plot shows 
the relative correction A^^ of Eq. (1291) due to the squark loops, summed over all quark 
flavors and with squarks degenerate in mass, and the lower plot the correction due to the 
gluino loop. The gluino loop correction exhibits a characteristic "dip" at the production 
threshold of a pair of gluinos, \^ = 2m^, reducing the LO cross section by about —10%. 
The corrections then increase with due to a ln(s/m?) dependence. A similar behavior is 
observed in the squark loop correction, which can reduce a^^ by up to about —4%, again 
due to a threshold at VI = 2mq, and then increases with ^/I. Both the gluino and the 
squark loop corrections do not depend on the squark mixing angles. 

Fig.[9](b) shows the influence of the vertex corrections of Fig.[3](b),(c), where in the upper 
plot it is assumed that m^^ = mi^ . The corrections are then independent of the top-squark 
mixing angle, because the term that is proportional to 9i cancels in Eq. (l23ll . The corrections 
show again a characteristic threshold behavior at VI = 2 nig. This threshold is introduced 
by the vertex correction containing the ggti,2 loop, shown in Fig. [271 This loop correction 
is the dominant contribution to the vertex corrections because of a large overall color factor 
of C\ = 9, arising due to the ggg coupling. If the top squarks are non-degenerate in mass, 
as it is the case in the lower plot of Fig. [9](b), the vertex correction depends very strongly 
on the stop mixing angle. Shown are corrections for m^-^ = 600 GeV and m^^ = 100 GeV 
for different values of the stop mixing angle. The largest correction occurs always at the 
gluino-pair threshold, VI = 2m^, for Of = tt/A and the smallest corrections are obtained 
for 6f = — 7r/4. For other mixing angles the corrections lie between these two curves. The 
corrections are furthermore the larger the more pronounced the mass splitting between m^^ 
and rrif^. The amount of allowed mass splitting, however, is restricted by the Ap parameter. 

The size of the box corrections is illustrated in Fig.fTOlfa) for m^-^ = 600 GeV and m^^ = 100 



19 



6 

4 
2 

-2 
-4 
-6 

f 7.5 
5 
2.5 


-2.5 
-5 
-7.5 
-10 
-12.5 



Gluon self energy: qq — > tt 



L m-fg-200GeV 
- m-'i g - 400 GeV 


1 1 1 1 1 1 1 1 _ 
Squark loops - 

- 


-400 600 800 


1000 1200 - 


J 1 — III — ^ 
^ m~ - 230 GeV 

— — H-WU VJC V 

: m?- 600 GeV 


^ — 1 1 1 1 ! 
Gluino loop 




1 , , , 1 , , , : 



Vertex Correction: qq — > tt 

' ' I ' ' ' ' I ' ' ' ' I ' ' 
m~ = 400 GeV 



8=; 




30 ^400 500 600 700 800 900 1000 1100 

25 

20 \- 

15 \- 

10 Y 

5 |-' 

- 

5 L = {600,100} GeV 



Sf = + Ji/4 

% = + 71/121 

9- = - Ji/4 i 




400 



600 



800 1000 
[GeV] 

(a) 



1200 1400 



350 400 450 500 550 600 650 700 750 800 
Vi" [GeV] 

(b) 



Figure 

vertex 



9: The relative correction Aqq{s) for qq annihilation due to (a) gluon self-energy and (h) 
corrections (upper plot: mi = = mi^). 



GeV and different values of the stop mixing angle. The upper plot shows the direct box 
contribution of Fig. [4](a) and the lower plot the crossed box contribution of Fig. [4](b). The 
threshold peak at \/I = 2mg originates from the two gluinos occurring in the Feynman 
diagrams of Fig. [H The corrections strongly depend on the stop mixing angle and in the 
threshold region are largest for 9^ = 7r/4, and can be suppressed very effectively for 9^ = 
— 7r/4. The direct and crossed box corrections have the same relative sign compared to 
the LO cross section, as discussed in Appendix [B] and are always negative in the threshold 
region. The direct box correction is usually three to four times larger than the crossed box 
contribution. Both box corrections depend on the masses of the spartners of the initial- 
state quarks. The corrections are small in this scenario, because of our choice of the squark 
masses, rriq = 2 TeV (for rriq > 1 TeV the corrections hardly change anymore). The direct 
box correction increases for smaller squark masses, e.g., for niq = 200 GeV at the gluino-pair 
threshold oi V§ = 460 GeV to Aqg = -22% (for rUg = 230 GeV, m^-^^ = 600, 100 GeV, % = 
+7r/4). However, the positive vertex correction also grows by the same amount and cancels 
the effect of the box corrections. We found that for every configuration with large box 
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corrections (and negative as always), the vertex corrections are positive and also large, 
so that both largely cancel, resulting into a small positive contribution to Agg. As already 



m 



i2' 



observed in the case of the vertex corrections, for top squarks degenerate in mass, m^^ 
the direct and crossed box corrections are independent of the top-squark mixing angle, and 
are largest for the smallest possible choice of their masses. 
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Figure 10: The relative correction Agq(s) for qq annihilation due to (a) box corrections and (b) 
when including the complete SQCD one-loop corrections for rUg = 230 GeV. 

One can now deduce from the detailed discussion above, for which MSSM scenarios the 
complete SQCD one-loop corrections, consisting of the sum of self-energy, vertex and box 
corrections (see Eq. ([4])), affect the qq — > ti total cross section the most. In Fig. [TOl(b). we 
show four representative scenarios, where we choose a light gluino mass of rrig = 230 GeV. 



If the top squarks are degenerate in mass, m^^ 



m 



t2' 



the corrections are independent of the 



squark mixing angle. If furthermore the stop masses are relatively light, i.e. smaller than 
500 GeV, no large corrections occur because of a cancellation between the positive vertex 
correction and the negative gluon self-energy. Only if the stops are very heavy, i.e. about 
2 TeV (dot-dashed line in Fig.fTOlfb)). the vertex and box corrections are strongly suppressed, 
and the negative gluon self-energy determines the NLO SQCD uu ti cross section. 
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If the stops are non-degenerate in mass, the corrections strongly depend on 6i (red solid 
line and blue dashed line in Fig. [TOlfb)). Large positive corrections can only occur, if the 
positive vertex correction becomes large (and all other, negative corrections are small) as 
is the case for 9i = n/A. Furthermore, the corrections are largest, if the lighter stop mass, 
mj^, is as light as allowed by the current mass limit, i.e. m^^ ~ 100 GeV, and for the 
largest possible mass splitting. The smallest correction for top squarks non-degenerate in 
mass occurs for 9i = — 7r/4 because then the vertex correction is negligible and the total 
correction is again dominated by the negative gluon self-energy. For all other values of the 
stop mixing angle, the relative correction Agg will lie between the red solid and blue dashed 
curves in Fig. [TOVb). 

2. Effects of NLO SQCD corrections in gluon fusion for unpolarized top quarks 

Top-pair production at the LHC is dominated by gluon fusion and, thus, the discussion 
of the SQCD one-loop corrections to gg — > ti at parton level is a good indicator of what to 
expect at the LHC. 

The impact of the gluon self-energy, top self-energy, and vertex corrections to the total 
partonic cross section in gluon fusion is shown in Fig. [TT], in form of the relative correction 
Agg of Eq. ( l29l ). The gluon self-energy correction to gg ti is displayed in Fig.fTTlfa). where 
the upper and lower plot shows respectively the squark-loop (summed over all quark flavors 
and with squarks degenerate in mass) and gluino-loop contributions. Both corrections are 
positive and much less pronounced than in qq annihilation. This is because the s channel 
contribution to the gluon fusion cross section (see Fig. [5]) reduces the total LO cross section 
by only 7%, interfering destructively with the t and u channels. For the same reason, the 
corrections are practically independent of the masses of the light squarks, since diagrams 
that contain u, d, c, s or b quarks occur only in the s channel correction. 

The upper plot of Fig. [TTlfb) shows the effects of the off-shell top quark self-energy 
contribution in the t and u channels of gluon fusion (see Figs. [5](b)). Again the SQCD 
one-loop corrections are the larger the lighter the gluino and the stops. They are positive, 
increase with rising up to about Agg = +3% at = 1200 GeV. They strongly depend 
on the stop mixing angle, if the squarks are non-degenerate in mass, and become largest 
for 6j: =11/4 (solid red line). A similar behavior can be observed, if m^^ = m^^ ^ii^d small, 
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Figure 11: The relative correction Agg(s) for gluon fusion due to (a) gluon self-energy, (b) (upper 
plot) off-shell top quark self-energy, and (b) (lower plot) vertex corrections. 

i.e. ~ 200 GeV. With increasing stop masses this correction quickly decreases. 

The effects of the vertex corrections are displayed in the lower plot of Fig. [TTlfb). The 
corrections are negative and exhibit the gluon-pair threshold behavior at = 2 due 
to the vertex diagrams including the ggii loop (see Fig. [27]). The largest corrections arise 
again for light gluino masses, a large mass splitting of the top squarks and small values of 
m^^. However, the influence of the mixing angle is small and changes the correction in the 
threshold region by maximally 1.2%. Again, if both top squarks become heavier the size of 
the corrections quickly decreases. 

The largest corrections in gluon fusion arise due to the box diagrams of Fig. [H as can 
be seen in Fig. [T2](a). In the gluino-pair threshold region the corrections can reach up to 



A 



gg 



-1-16% for light gluino masses of 230 GeV, large mass splitting of the top squarks 
and 9t = 7r/4. The corrections depend again strongly on the stop mixing angle, if the stop 
masses are not equal, as displayed in the lower plot of Fig. [T2](a). If m^^ = m^-^, the box 
corrections at the gluino-pair threshold reach Agg = +8% for stop masses of 200 GeV and 
only Agg = +2% for masses of 500 GeV. 
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Since the top self-energy and vertex correction add destructively and the gluon self-energy 
is negligible, the complete SQCD one-loop correction to gluon fusion is dominated by the 
box corrections and exhibits the same characteristics. This is illustrated in Fig.[T2](b). where 
we show the impact of the complete SQCD one-loop correction to the partonic total cross 
section of gluon fusion for four different MSSM scenarios, which differ by the choices of top- 
squark masses and mixing angles. The corrections are practically independent of the light 
squark masses, m^, because they contribute only to the gluon self-energy. The corrections 
are mostly positive, with a size of maximally = +15% at the gluino-pair threshold. 
They can become negative, reaching maximally A^^ = —2.5% for ~ 1 TeV. 
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Figure 12: The relative correction Agg{s) for gluon fusion due to (a) box corrections (upper plot: 



mi 



mr 



mi ) and (b) when including the complete SQCD one-loop corrections. 



3. Hadronic cross sections to unpolarized pp,pp tt at NLO SQCD 

The effects of the NLO SQCD corrections to the hadronic cross sections result from the 
combination of the already discussed effects to the partonic cross sections and the PDFs (see 
Eq. dl])). The hadronic ti cross sections at the Tevatron Run II and the LHC are dominated 
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respectively by the quark and gluon PDFs, which both emphasize parton-level effects in the 
vicinity of the ti threshold, = 2mj = 350 GeV, but are rapidly decreasing for increasing 
values of VI. The choice of the values for the MSSM input parameters, rrig, rrif.^, mi^ and 9i, 
is guided by the discussion of the SQCD effects at the parton level of Sections IlII A 11 [Til A 2[ 
We first study the dependence of the relative correction A of Eq. ( 1301) on the gluino and 
heavier stop masses for four different choices of the light stop mass and the stop mixing angle. 
These choices are representative for the possible choices that affect the ti cross sections the 
most. 

The impact of the SQCD one-loop corrections on the total hadronic cross section in 
dependence of the gluino mass is shown in Fig.fTSVa) for the Tevatron Run II and in Fig.[T3l(b) 
for the LHC. The corrections are largest for gluino masses of 200 — 250 GeV, since for these 
masses the gluino-pair threshold lies in the vicinity of the ti threshold. For rrig = 230 GeV 
the SQCD one-loop corrections vary between A = —4% and A = +5% at the Tevatron Run 
II, and between A = —1.5% and A = +5.5% at the LHC. They decrease to approximately 
A = —0.5% (Tevatron Run II) and A = +1% (LHC) for gluino masses of 500 GeV. Since 
large corrections always originate from the gluino-pair threshold, the relative corrections 
decrease for heavier gluinos. 



pj5 tt (VS = 1960 GeV) pp tt (Vs = 14TeV) 




m; [GeV] m- [GeV] 

(a) (b) 
Figure 13: The relative correction A due to SQCD one-loop corrections in dependence of the gluino 
mass, rrig, at (a) the Tevatron Run II and (b) the LHC. 

The influence of the top-squark masses on the SQCD one-loop corrections is shown in 
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Fig. [T4l(a) for the Tevatron Run II and in Fig. [T4l(b) for the LHC. Varied for these plots is 
either the heavier top-squark mass alone or, in case of the dot-dashed line, both stop 

masses are varied together. The solid red, dashed blue and dotted black lines show that 
the corrections increase with increasing difference in the top-squark masses, as discussed 
at the parton level in Section IIII A 11 IIII A 2[ The solid red and dashed blue lines end at 



mi 



650 GeV and the dotted line at m^^ ^ 850 GeV, because for larger values of m^^, the 



7l|. At the Tevatron Run II, if the top 



parameter Ap exceeds the current limit of 0.0035 
squarks are degenerate in mass, the largest correction arises for very large top-squark masses 
because then the correction is determined by the gluino loop in the gluon self-energy, which 
amounts to A ^ —4% at hadron level. For > 2 TeV the correction does not increase 
much anymore and stays at about A = —4%. At the LHC, the gluon self-energy plays a 
very minor rule, and all other corrections depend on the top squark mass. Therefore, the 
corrections decouple and approach zero if both top squarks are heavy. 
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Figure 14: The relative correction A due to SQCD one-loop corrections in dependence of the heavier 
top-squark mass, mi_^, at (a) the Tevatron Run II and (b) the LHC. 



The relative corrections to the and pt distributions of Eq. ([30l) are shown in Fig. [15] 
and Fig. [TBI respectively. Both the Mu and the px distribution exhibit the character- 
istic effects of the gluino-pair threshold, which can lead to a significant distortion of 
the shape of these distributions. At the Tevatron Run II, for instance, when choosing 
"^ti,t2 ~ 600' -'-00 GeV, 9t = n/i and nig = 230 GeV, the corrections to the Mtiipr) distri- 
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Figure 15: The relative correction A(M(f) due to SQCD one-loop corrections at (a) the Tevatron 
Run II and (b) the LHC. 

bution increase for M^j < 460 GeV(pT < 150 GeV) and can enhance the LO distribution by 
up to A{Mtt){A{pT)) = +15(+9.5)% and then decrease quickly for larger values oi Mttipr)- 
This behavior suggests that the relative corrections to the total tt cross section can be en- 
hanced by applying cuts on the px of the top quark. For instance, when restricting the top 
quark px to the range 75 GeV < pr < 170 GeV, the relative correction at the Tevatron 
Run II reaches A = +7.1% (instead of A = +5.4% without cuts, see Fig. [13]). Keeping the 
same top-squark parameters but choosing rrig = 260 GeV and 100 GeV < pr < 210 GeV, 
we find A = +5.5% (instead of A = +4.1% without cuts). At the LHC, choosing 
mi^^l^ = 600, 100 GeV, 9^ = 7r/4, nig = 230 GeV and 100 GeV < Pr < 170 GeV, we 
find A = +7.5% (instead of A = +5.7% without cuts). We also studied the effects of rapid- 
ity cuts and found that they do not affect the relative corrections much, because the NLO 
SQCD corrections are quite stable with respect to the rapidity distribution of the top quark. 

In view of anticipated experimental uncertainties of about 10% and 5% in the total ti 
production rate at the Tevatron Run II and the LHC, respectively, and a current theoretical 
uncertainty of the QCD prediction of about 12% Q, [ll|, it will be difficult to reach sensitiv- 
ity on SQCD effects in this observable. The and px distributions exhibit an interesting 
signature of SQCD one-loop corrections, but they are strongly affected in only a few bins 
and the final verdict on its observability must be left to an analysis, that includes for in- 
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stance top decays and the detector response. However, the inclusion of the known SUSY 



electroweak one-loop effects 



47| may help to enhance SUSY loop-induced effects. First re- 



sults of combined SUSY electroweak and SQCD one-loop corrections have been presented 
at the Tevatron in Ref. [s^, but a detailed study at the LHC and the Tevatron Run II still 
needs to be done. 
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Figure 16: The relative correction A(j)t) due to SQCD one-loop corrections at (a) the Tevatron 
Run II and (b) the LHC. 
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B. Polarized top-pair production at NLO SQCD 

In this section, we consider the production of polarized top-quark pairs at the Tevatron 
Run II and the LHC. We are working in the helicity basis of left- and right-handed top 
quarks, as discussed in Section III A[ The incoming quarks and gluons are still considered 
to be unpolarized. In the following, we use the notation that the indices 'L' and 'R' denote 
left {\t{t) = —1/2) and right-handed {Xt{t) = +1/2) top(antitop) quarks, respectively. 

Since the SQCD couplings of gluinos and squarks contain an axial-vector part (see Ap- 
pendix [Al Fig. [24j), the SQCD one-loop corrections may affect the production of left and 
right-handed top quarks differently. To study these differences in detail, we first discuss in 
Sections IlII B II {qq annihilation) and Section [ill B 21 (gluon fusion) the impact of the SQCD 
one-loop corrections on the total partonic ti cross sections of Eq. ^ for each top and antitop 
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helicity state separately, using the following relative corrections at the parton level: 



Xab _ ^qq,99 (^*' ~ ^gg-gg(^*' 

(^qq,99[M, -\t) 



with ab = LL, RR denoting the relative correction for spin-like ti production with At = 
— 1/2, Af = —1/2 and Aj = +1/2, Aj = +1/2, respectively. Similarly, ab = LR,RL denotes 
the relative correction for spin-unlike ti production with \t = — l/2,Af = +1/2 and \t = 
+1/2, Af= —1/2, respectively. Naturally, since QCD preserves parity, there is no difference 
in the spin-unlike amplitudes at LO QCD, i.e. al^Pgg{+l/2, -1/2) = a^9gg{-l/2,+l/2). 
Since we are not considering any CP violating couplings, the spin-like amplitudes are the 
same at both LO QCD and NLO SQCD, i.e. a^^^f ^'^(-1/2, -1/2) = a^^^f ^'^(+1/2, +1/2). 

Differences in the spin-unlike amplitudes at NLO SQCD manifest themselves in parity- 
violating polarization asymmetries, which we study in detail at hadron level in Section llll B 31 



12|: 



We consider the following differential and integrated polarization asymmetries 

• The left-right asymmetry in the Mti distribution 

dALR{Mtt) = 1 jj—- (32) 

daRL/dMu + daLRl dMu 

and in the total hadronic ti cross section 
. ctrl-ctlr , . 

^LR = ■ , [OO) 

CTRL + CTLR 

where we introduced the notation daiR^RL) = da^Loi^t = '~l/2(+l/2), Af = 
+l/2(-l/2)) with the hadronic NLO SQCD cross sections of Eq. ((3]). 

• The left-right asymmetry in the distribution, when assuming that the polarization 
of the antitop quark is not measured in the experiment (denoted by f/=unpolarized), 



SAiMu) 



{daRL/dMu + daRR/dMu) - {daLL/dMu + daLRl dMti) 
{daRL/dMu + daRR/dMu) + {daLL/dMa + daLR/dMu) 
daRu/dMu - daLu/dMu 



daNLo/dMu 

and the corresponding asymmetry in the total hadronic ti cross section 



(34) 



^ _ (o'flL + cfrr) - {<^LL + (^lr) _ cfru - am ^^^^ 

(^RL + (^RR + (^LL + <7lR (^NLO 

with the unpolarized NLO SQCD cross sections, da^LO = Tl,\t At=±i/2 da]siLo{S, Xt, Af). 
In this way, one complicated polarization measurement is avoided. 
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These parity-violating asymmetries are zero for stop mixing ang les 61= ±\ (then A/ = 0, 
see Appendix [A]) and top squarks degenerate in mass. The latter is due a cancellation of 
the parity- violating terms in the sum of the ti and ^2 contributions, since = — Aj and for 
equal stop masses the terms multiplying A-^ are the same. 

In addition to the parity- violating polarization asymmetries, we study in Section IIII B 31 
the spin correlation functions of Ref. [76|, which are parity-conserving asymmetries in the 
spin-like and spin-unlike contributions to the ti cross sections. They are defined as follows: 
^ _ {daRR/dMu + daLL/dMu) - {daLn/dMu + daRL/dMa) 



{daRR/dMu + daLL/dMtt) + {daLR/dMu + daRi/dMu) 
{daRR/dMtt + daLL/dMu) - {daLn/dMu + daRL/dMa) 



daNLo/dMu 

and 



(36) 



^ _ {(^RR + CTll) - (O'LR + Q-Rl) _ {(T RR + CTll) " (cTlR + CT Rl) ^^^^ 
{(^RR + (^Ll) + {(^LR + (^Rl) (^NLO 

1. Effects of NLO SQCD corrections in qq annihilation for polarized top quarks 

In this section, we study the effects of the SQCD one-loop corrections on the partonic 
cross sections for polarized top quarks in qq annihilation. We are especially interested in the 
differences between the corrections to t^t/j and tRtL production, since they will determine the 
size of the polarization asymmetries at hadron level at the Tevatron. As in the unpolarized 
case, we only show results for uu ti, representatively for all gg-initiated processes. In 
Fig. [I7l we show the relative correction of Eq. (l3Tl) for the sum of the direct and crossed 
box corrections, separately for each top(antitop) polarization state. The characteristics 
of the box corrections are similar to the ones observed in the unpolarized case, which we 
discussed in Section IIII A 1[ In particular, the polarized box correction also exhibits the 
threshold behavior at = 2mg. Fig. [TWa) shows the correction in dependence of the 
partonic CMS energy for a small gluino mass of nig = 230 GeV and for a choice of the stop 
mixing angle, where we expect the left-right asymmetry to be large. The dependence of the 
box correction on the stop mixing angle is shown in Fig. flTlfb) with chosen to be at the 
gluino-pair threshold. As can be seen, the differences between the relative corrections to 
t^tR and t/jti production are largest for 6^ = ±7r/2, 0. 

In Fig. [18] we show the relative correction due to the complete SQCD one-loop correction 
in dependence of the partonic CMS energy (upper plots) and the stop mixing angle (lower 
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Figure 17: The relative corrections AJ^g(s) {ab = LR, RL, RR = LL) due to the box correction for 
polarized top-quark pairs in qq annihilation in dependence of (a) the partonic CMS energy and (b) 
the stop mixing angle. 

plots). Since the box correction is very small compared to the complete correction, and 
the gluon self-energy is independent of the top(antitop) polarization (and thus cancel in 
the polarization asymmetries), the differences between the relative corrections to tLtn and 
tntL production observed in Fig. [18] are mainly due to the vertex correction. Fig. [TSVa) 
shows the corrections for rrig = 230 GeV and m^-^ = 100 GeV, and Fig. [TSlfb) the corrections 
for somewhat heavier particles, rrig = 300 GeV and m^^ = 250 GeV. The vertex correction 
introduces large partonic left-right asymmetries in the gluino-pair threshold region and at 
high for 6f = ±7r/2,0. Furthermore, as illustrated in Fig. [18] (lower plots), the largest 
relative corrections occur in spin-like top-pair production (A^^ = A^^) at 9^ = ±j. In this 
case, the relative difference between the corrections to spin-like {tj^R^tLti) and spin-unlike 
{tRiL,ti,tfi) top-pair production increases with increasing gluino mass. We found that if the 
stops are degenerate in mass, the corrections to spin-like and spin-unlike ti production are 
of comparable size. If in addition the top squarks are very heavy (m^^ = m^^ ~ 2 TeV), the 
corrections are completely independent of the top polarization states, because then they are 
determined by the gluon self-energy, as discussed in Section [III A 1[ 

2. Effects of NLO SQCD corrections in gluon fusion for polarized top quarks 



When studying the effects of the SQCD one-loop corrections on the partonic cross sections 
for polarized top quarks in gluon fusion, we again are especially interested in the differences 
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Figure 18: The relative corrections Aq|(s) (a6 = LR, RL, RR = LL) due to the complete SQCD 
one-loop corrections for polarized top-quark pairs in qq annihilation in dependence of the partonic 
CMS energy (upper plots) and the stop mixing angle (lower plots) for two sets of choices of the 
MSSM parameters: (a) m-g = 230 GeV, m^^ ^^ = 600,100 GeV and (b) m-g = 300 GeV, m^^ ^^ = 
800,250 GeV. 



between tf>iL and tiiji production, since this will determine the size of the polarization 
asymmetries at hadron level at the LHC. As for qq annihilation, also in gluon fusion the gluon 
self-energy correction is independent of the top(antitop) polarization states. In Fig. [TWa) 
we therefore show the relative corrections A^^ of Eq. (l3Til separately only for the top self- 
energy and vertex corrections, choosing nig = 230 GeV and m^-^ ^ = 600, 100 GeV. As can be 
seen, the top self-energy and the vertex corrections depend strongly on the polarizations of 
the top and antitop quarks, inducing differences between the corrections to tjitL and tii^ 
production of several percent. However, the partonic left-right asymmetries due to these 
two corrections are of opposite sign and therefore almost cancel in the complete SQCD 
one-loop correction. Also the corrections to spin-like top-pair production have opposite 
signs and interfere destructively. Therefore, the relative correction due to the complete 
SQCD one-loop correction, shown in Fig.fTWb). is largely determined by the box correction. 
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Figure 19: The relative corrections A^^(s) {ab = LR, RL, RR = LL) for polarized top-quark pairs 
in gluon fusion due to (a) top self-energy and vertex corrections and (b) the complete SQCD one- 
loop correction in dependence of the partonic CMS energy. Also shown is the stop mixing angle 
dependence of the complete SQCD one-loop correction ((b) lower plot). 

The relative correction for spin-like top-pair production is much larger than the one for 
tRtLytLtR production, which is still the case for larger values of rrig. The difference in the 
box correction to tutL and titi? production is relatively small and therefore also the partonic 
left-right asymmetry induced by the complete SQCD one-loop corrections is small. 



3. Hadronic cross sections to polarized pp,pp — > tt at NLO SQCD 

In the following discussion of parity violating effects in polarized ti production at the 
Tevatron Run II and the LHC, we benefit from the detailed study of these effects at the 
parton level in Sections IIII B 11 IIII B 21 so that we can restrict our discussion to the following 
two sets of choices for the stop masses and mixing angle: 

(/) = 600 GeV, nii^ = 100 GeV, % = n/U 
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{II) mi^ = 800 GeV, m^-^ = 250 GeV, % = (38) 

and to a few values of the gluino mass. These parameter sets are representative for the 
choices that yield the largest numerical impact. The top-squark mixing angle of set (I) is 
chosen to be 7r/12 instead of zero, since in this case smaller values of 9^ render the soft 
supersymmetry breaking mass M?^ negative. 

At the Tevatron Run II, the integrated left-right asymmetries in the total hadronic cross 
section of Eqs. ([33]), ([35]) are largest for set (I): We find Alr = -1.08%,^ = -0.79% for 
m~g = 230 GeV, and Alr = -0.72%, A = -0.53% for nig = 300 GeV. These small asym- 
metries are clearly not observable at the Tevatron: For instance, the statistical significance 



Ns of Eq. (19) in Ref. [12| only amounts to Ng = 2.0(1.7) for \Alr\ = 1.08%(|^| = 0.79%) 
for an integrated luminosity of £ = 8fb~^, while conservatively Ns > 4 is required to 
be statistically significant. In Fig. [201(a) we show the parity- violating asymmetries in the 
Mtt distribution, SALR{Mtt) of Eq. §2) and SA{Mtt) of Eq. ([34]), for sets (I) and (II) and 
a gluino mass of 230 GeV and 300 GeV. The asymmetry 6A{Mtt) could be an interest- 
ing observable at the Tevatron, since the cross sections for the spin-like polarization states 
{tLti, iRtR) are much smaller than the ones for the spin-unlike polarization states (tL^R, ^r^l) 
and, thus, do not significantly decrease this asymmetry compared to 5ALR{Mtt). Due to 
characteristic peaks at the gluino-pair threshold, 5A{Mtt) can reach — 3.4%(— 4.1%) for set 
(I) and rrig = 230(300) GeV, and — 2.1%(— 3.0%) for the heavier top-squark mass set (II). 
If the polarization of the antitop quark is also measured, the resulting left-right asymmetry 
6Alr{Mu) of Eq. ([32]) could reach -4.7% (set (I) and m-g = 230 GeV) and -4.6% (set (I) 
and mg = 300 GeV). 

In Fig. [201(b). we show the parity- violating asymmetries in the Mu distribution of 
Eqs. ( [321) . ( l34l ) at the LHC. Since at the LHC the cross section oi tRtR^titL production is 
much larger (roughly three times) than the one for tiiji^tjiiL production, the asymmetries 
6A{Mtt), Aof Eqs. are significantly smaller than SAlr{Mu),Alr of Eqs. ([32]). ([33D. 

At the LHC, the left-right asymmetry 6ALR{Mtt) reaches a maximum for set (I) of —2.0% 
for nig = 230 GeV and —1.1% for rrig = 400 GeV. In comparison, if the polarization of the 
antitop quark is not measured, we find for the same set of MSSM parameters that 6A{Mtt) 
is maximally —0.64% and —0.58%. The asymmetries in the total hadronic cross section, 
Alr. of Eq. f[33l) and A of Eq. f[35ll . reach for set (I) a maximum value of —0.24% and 
-0.08% for m-g = 230 GeV and -0.37% and -0.13% for m-g = 280 GeV. For set (II) the 
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Figure 20: Parity-violating asymmetries in the M^i distribution for polarized top-quark pairs 
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plots) and for polarized top and unpolarized antitop quarks (lower plots) at (a) the Tevatron Run II 
and (h) the LHC. 

asymmetries are maximal at roughly rrig = 300 GeV with Alr = —0.26% and A = —0.09%. 
For both parameter sets, Alr and A decrease slowly for heavier gluinos. At the LHC, 
integrated parity-violating asymmetries as small as 0.1% may be accessible: For instance, 
we find statistical significances of Ns = 9.0(5.3) for \Alr\ = 0.37%(|^| = 0.13%) (with 
£ = 30fb-^). 

In the remaining part of this section we will provide numerical results for the spin cor- 
relation functions of Eqs. ( l36l) . (i37l ) at the LHC. Fig. [2T](a) shows the difference between 
the LO QCD and NLO SQCD predictions for the spin correlation function C of Eq. fl36l) 
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in dependence of the top-squark mixing angle. A comparison of the LO QCD and NLO 
SQCD predictions for the spin correlation function C of Eq. (l37l) is shown in Fig. [2T](b). 
The LO QCD predictions are obtained from Eqs. ( l36ll . (l37l) by replacing the NLO cross sec- 
tion with the LO cross sections. As shown in Fig. [TOlfb). the SQCD one-loop corrections to 
the tRin^tLtL and tLtR.tRtL production can differ considerably, resulting in large deviations 
from the LO QCD prediction for C, and thus can lead to potentially observable effects at 
the LHC. Since this difference strongly depends on the stop mixing angle, this spin corre- 
lation function may have the potential for extracting information about 6i. In Fig. [22l we 




% [rad] M„ [GeV] 
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Figure 21: Spin correlation functions in polarized ti production at the LHC at LO QCD and NLO 
SQCD. Shown is (a) the difference between the LO QCD and NLO SQCD spin correlation functions, 
Co andC, in dependence of the top-squark mixing angle % and (b) C in dependence of the invariant 
ti mass. 

discuss the difference between the LO QCD (Co) and NLO SQCD (C) predictions for the 
spin correlation function of Eq. (1371) in more detail. As can be seen in Fig. [221(a). a deviation 
from the LO QCD prediction of maximal 3.7% can be achieved for light stop and gluino 
masses and a top-squark mixing angle of 9f = n/A. Since the spin correlation functions 
C, C do not suffer from luminosity uncertainties (they cancel in the cross section ratios), 
and due to the high tt yield at the LHC, they may be interesting observables to search for 
loop-induced SUSY effects in ti production. To see how these effects compare to theoretical 
uncertainties in the predictions for C we show again in Fig. [221(b) the difference Cx — Cq 
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with C X = C for m^^ = 600 GeV, m^^ = 100 GeV, 9^ = n/i and a light and heavier gluino, 
together with the uncertainty bands induced by the LO QCD scale dependence, varied be- 
tween 1711/2 < jip = jJ'R < 21711, and an experimental top-mass uncertainty of Arut = 1 
GeV. The bands are calculated with Cx being the LO spin correlation function for different 
values of rut and ^r, ftp, and Co is calculated ioi rrit = fip = I^r. = 175 GeV. Since the scale 
uncertainty is obtained at LO QCD, we expect the corresponding band to be considerably 
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less pronounced when including the known QCD NLO results 
view of an anticipated relative experimental error on C of about 4% 
has the potential to be sensitive to loop-induced SUSY effects. 
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In 



77||, this observable 
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Figure 22: Difference between the LO QCD and NLO SQCD spin correlation functions, Co and C, 
at the LHC in dependence of the invariant tt mass. In (b) we compare the NLO SQCD prediction 
for this difference (dotted lines with Cx = C and rrif:^ = 600 GeV,mj:^ = 100 GeV,9^ = 7r/4, = 
230,400 GeV) to the LO QCD prediction with Cx describing the LO spin correlation function for 
varying top mass and renormalization and factorization scales. The bands are obtained due to a 
variation ofrtit by ±1 GeV (grey, cross-hatched band with the lower line corresponding to rrit — l GeV 
and the upper line to rrit + l GeV) and due to a variation of fip and fiR by mt/2 < fip = i^ir < 2mt 
(red, hatched band). Cq in (a) and (b) is calculated for mt = hf = IJ'R = 175 GeV. 
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C. Comparison with existing calculations 



The NLO SQCD corrections to qq tt for unpolarized top quarks have also been cal- 



culated in Refs. 



53 



54 l55|, however in Ref. [53|| the contribution of the gluon self-energy 



and the crossed box diagram (Fig. [4](b)) is mis sing . We compared our results analytically 
and numerically with those presented in Refs. [sd IssI and found agreement, if we adjust 
for a missing (—1) sign for the direct box (Fig. [4](a)) contribution in Ref. [5J| and missing 
(— 1) signs for the direct (Fig. [4](a)) and crossed box (Fig.[4](b)) contributions of Ref. [s^. A 
detailed discussion of the box contributions can be found in Appendix [B1 The NLO SQCD 
corrections to gluon fusion for unpolarized top quarks have also been calculated in Ref. [SOj. 
When we use their choices of the SM and MSSM input parameters and assume a. = 0.120 



we were able to reproduce the numerical results presented in the figures of Ref. 



50|. 



IV. CONCLUSION 



We studied in detail the effects of SQCD one-loop corrections to the main production 
processes for polarized and unpolarized strong tt production, qq tt and gg tt, at the 
Tevatron Run II and the LHC. We presented numerical results for the total tt production 
rate, the invariant tt mass and top transverse momentum distributions in unpolarized ti 
production and for a number of asymmetries in polarized tt production for different choices 
of the MSSM input parameters, m^, m^^, m^^ and 6i. We found that the largest corrections 
occur for the lightest gluino allowed by current experimental limits and large stop-mass 
splittings. For instance, for m^-^, m^^ = 600,100 GeV and 6i = 7r/4, the SQCD one-loop 
corrections enhance the LO total hadronic ti cross section by 7.1% (when restricting the 
top quark px to 75 GeV < pr < 170 GeV) at the Tevatron Run II for nig = 230 GeV and 
by 7.5% (100 GeV < pr < 210 GeV) at the LHC for m-g = 260 GeV. The NLO SQCD Ma 
and Pt distributions can be significantly distorted due to a gluino-pair production threshold 
at = 2mg. For and pt values in the vicinity of this threshold, the SQCD one-loop 
corrections can enhance the LO and pt distributions by 15% and 9.5%, respectively. 
When considering polarized ti production, we studied the effects on parity-violating and 
parity-conserving asymmetries in the total hadronic cross section and the distribution. 
We found that in view of the anticipated top-quark yield at the Tevatron Run II, it will 
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be difficult to observe loop-induced SUSY effects in polarization asymmetries. At the LHC, 
however, we find promising effects: The parity-violating asymmetries in the production of 
left and right-handed top and antitop quarks can reach up to \6Alr{Mu)\ = 2.0% in the Mu 
distribution and \Alr\ = 0.37% in the total hadronic cross section. The parity-conserving 
spin correlation function, C, that describes an asymmetry in the Mu distributions of spin- 
like and spin-unlike ti production can differ from the LO QCD prediction by up to 3.7%. 
These effects are promising enough to motivate a future study which includes top-quark 
decays and the detector response, in order to determine whether they are observable at the 
LHC. Such a study should also include SM and SUSY electroweak one-loop corrections, 
since they can either enhance or diminish the SQCD induced effects. 
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Appendix A: FEYNMAN RULES 

The relevant QCD Lagrangian of the MSSM is given in Ref. with the convention of the 
covariant derivative -D^ = 9^ + igsG'tT"- with the gluon field G^, strong coupling parameter 
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Qs = i/ivraj, and T° = A°/2, where A" denote the Gell-Mann-matrices. Furthermore the 



methods of Ref. 
In 



66| are apphed to address Feynman graphs containing Majorana particles. 



66l | a "fermion flow", denoted by a thin line with arrow, is introduced in addition to the 
standard "fermion number flow" to deal with diagrams containing fermion number violating 
fermion chains. The relevant SQCD Feynman rules used in this paper are given in Figs. \23l 



~*k 
1i. ^ 



^9sT^l{p-p)^,5^J 




1i 



*k 



,2fSa 



ij 



Figure 23: The Feynman rules for the triple and quartic squark-gluon interactions. a,b,c = 1 ... 8 
and k,l = 1...3 are color indices and i,j = L,R and i,j = 1,2{L,R) with(without) mixing are 
the squark indices. Here d"'''' is total symmetric, defined by d"^" = 2Tr [{T" ,T^}T'']. 





{gi - gpi^) 



Figure 24: The g-qj-q - vertex written in generalized form in terms of scalar and pseudo scalar 

L R Li Ft L Ft 

couplings, QsJo , with gs' = ±1 and gp' =1. The thin line with arrow represents the fermion 



flow of Ref. 



66| 



When considering g^-g/j-mixing the interaction eigenstates are replaced by the mass eigen- 
states in the interaction Lagrangian, as discussed in Section III CI which has the following 
impact on the squark couplings of Fig. | 



g'/ = cose,-^?f"^±sine,-r7f'^ 
g'/ = cos^,-^^'^±sin^,-^?f'^ 



For convenience we introduce the abbreviations 



\f = (gi)' ± (glY and \f = 2gigj, . 



(Al) 



(A2) 



The SM QCD Feynman rules which we need to calculate the tt cross sections are given in 

Fig.n 
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Figure 25: The g-g-g - vertex with the SU(3) structure constants fabc, and the gluino and squark 
propagators. 




iS, 
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9sfaia2aa X 

{(fci - k2)^,g^,f,2 

+ {k2 - fc3)Ml3M2M3 

+ (^3 - kl)^l.29^^3^l^l} 



fi,a 



P 



9, 



gsPf,f 



abc 



Figure 26: SM QCD Feynman rules. The long thin lines with arrows denote the fermion flow of 



Ref. 



66| 



To obtain the Feynman rules with reversed "fermion flow" for diagrams containing Dirac 
fermions (Fig. [241 gqq-vertex and quark propagator in Fig. [26l) one has to replace in the 
above rules the strings of Dirac matrices Fj = 1, ^75, 7/^75, 7^*, cr^^ by 



T[ = r]iT, = CTjC-' with 



(A3) 



1 for Ti = 1, 275, 7^75 
-1 for Fi = 7^, a^^ 
where C denotes the charge conjugation operator. Further rules for dealing with external 



spinors can be found in Ref. 



66|. 



Appendix B: ANALYTICAL NLO SQCD CORRECTIONS TO qq tt 

With the Feynman rules presented in Appendix A and the counterterms of the renor- 
malization procedure as described in Section III Bl the SQCD one-loop corrections to the 
qq annihilation subprocess can be given in the compact form of Eq. (Jll). The explicit ex- 
pressions of the UV flnite (after renormalization) gluon self-energy contribution, n(s), is 
given in Eq. (fT2l) of Section III B[ The renormalized gqq vertex is paramatrized in terms of 



UV finite (after renormalization) formfactors, Fy, Fm, Ga (introduced in Ref. [29|), describ 



ing the vector, magnetic and axial vector parts, respectively. The SQCD one-loop vertex 
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corrections of Fig. [27] modify these formfactors as follows: 



a 
An 



—Fv(s,m^ 



.3=1,2 



as 
Air 



Fuis^mq) 



- Xj2mgmqCo{s,mg^,mg,mg)]} + 6Zv 



3 

47r 



(Bl) 



- m-g mg{Co + Ci + C2) (s, m^, m,^. , m^^. ) 



2 ^ 



777 - 

XUCi + C2 + Cn + C22 +2Cu) + K — {Ci + Q 



m„ 



(B2) 

[s,mg^,mg,mg 



—GA{s,mg 
47r 



Q: f 1 3 

^ \ Q^fCoo{s,mg,mg^,mq^) + -Xf [-Bo{s,mg,mg 



+ (2Coo -(m^ -^g + ^i)Co - 2m^(Ci + C2))(s, mg-^., m^, m^)]! + 6Za (B3) 

with 5Zv,A of Eq. ffTSl) and A^, A^ defined in Eq. (IA2I) . The two point functions, 

Bo{s,mi,m2) = Bo{s,ml,ml), and the three point functions, [Ci,Cim]{s,mi,m2,m3) = 
[Ci, Cim]{m'^, s, m^, mf, m^, m^), follow the notation of Ref. jfi^. 






Figure 27: Renorm,alized gluon- quark- quark vertex correction at NLO Susy-QCD. Graphs containing 
squarks are summed over the squark mass eigenstates j=L,R(l,2) (no(with) mixing). 



I 9 



t 




(a) (b) 
Figure 28: (a) Direct box diagram and (b) crossed box diagram contributing to qq — > ti at NLO 



SQCD, where we explicitly indicate our choice for the fermion flow (see Ref. f6d]) 



Finally, we provide the explicit expressions for the direct (Bt) and crossed box diagrams 
(Bu) in qq annihilation. To determine the relative sign between the direct and crossed box 
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diagrams (as well as between box and Born contributions) we apply the rules of Ref. [66|. 
Fixing the reference order as tiqq and choosing the fermion flow as indicated in Fig. [28] 
we only need to assign an additional minus sign to the direct box contribution. Because of 
differences in previous calculations of the direct and crossed box diagrams in qq annihilation, 
we also provide the explicit expressions for the corresponding matrix elements, SA4lg^ and 
^■^box^ of the Born matrix element 

^MZx ) ^"^M ii-^-^3+mg) {Tl'l^) 'vr {P2)un (pi ) ^Fg^^ J) i 4 - ^ + m^) F2'j^ M„ (^4) 

Ml'= Ur{P2)H9sVnY)Vn{pi) {-l^-f5''')MP^){-WsTLl'')uM (B4) 

withX(pi) = (27r)''[A;2-m|][(p4-A;)^-m?][(p3 + A;)2-m?][(pi + A;-p4)2-m|], s = {p^+Pif 
and z, j = 1, 2(L, R) denoting squark indices, and a, b, k-r color indices. The couplings ^'^{^^ 
are specifled in Fig. [24] and the momenta of the external particles are defined in Fig. [H The 
matrix element for the direct box diagram, 6M.l^^, of Eq. (IB4I) receives an additional minus 
sign because of the permutation parity of the spinors with respect to the Born diagram of 
Fig. [T] that has been chosen as the reference diagram. 
Direct box 

The Feynman diagram of the direct box is shown in Fig. [28](a). The corresponding spin-like 
contribution to the NLO SQCD matrix element squared of Eq. ((41) reads (with z = cos 9) 

Bt{s, i Xt = 1/2, Xt = 1/2) = Bt{s, i Xt = -1/2, A* = -1/2) = ^ !^(l - 

X I Af'^mt [2m^tiDo + 20^ + D22 + 2A3) - 2sl3^^{Di + D^^) - S(l + A')/^ii - 4Doo] 

+ A'^m-g [4m^(Do + D2) - 2sA'Di] + Af''2mlmtDo^ (B5) 

with [Di, Dim] = [Di, Dim]{mf, 0, 0, m^, i, s, m? , m?, m? , m?) defined in Ref. 6^ and depend- 
ing on the squark fiavors i and j. The spin-unlike parts of the direct box can be written as 
follows 

Bt{s,i, Xt = 1/2, Xt = -1/2) = {Bt + B',){s,i) 

Bt{s,i,Xt = -1/2, Xt = 1/2) = [B^^ - B\){s,i) (B6) 
with 



-ml(3tAzD^ + ^m?(l - 2(3tz + z''){2D2 + Dq + D22) 
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Bl{s,i) = [-2m?/3,(l + ^2)(Di + Di2)-(A + 2^ + A^')K'^ii + ^oo) 



16 



+ A'^Al^m-gmt [-2^(1 + z^)D^ - (A - 2z + Az2)(Z}2 + D^)] 
+ ^Af^m|(A-2^ + A^2)^^| 

The coefficients are defined by tlie mixing matrices of tlie squarks occuring inside tfie 
loop. In case of a real mixing matrix these coefficients can be expressed in terms of the 
squark mixing angles 9q and % : 

1 cos^ 6q cos^ 6't + sin^ Oq sin^ if j 
A\ = {-iysm29i, 
Af = - (-l)' cos2^g , 



^3 



COS 6'g cos 9^ + sin 6*, sin 9^ if z = j 
cos^ 9q sin^ 6'j + sin^ 9q cos^ 6*^- if i 7^ j 

Af' = - [(-ly cos2%+ (-l)^cos2eg] 

Af'' = - [(-ly cos2^,-- (-l)^cos2^,-] (B8) 
Crossed box 

The Feynman diagram of the crossed box is shown in Fig. [281(b) . The corresponding spin-like 
contributions, Bu{s,i, Xt = 1/2, Af= 1/2) = Bu{s,i, \t = — l/2,Af = —1/2), and spin-unlike 
contributions, Bu{s,i,Xt = 1/2, Xt = -1/2) and Bu{s,i,Xt = -1/2, Xt = 1/2), to the NLO 
SQCD matrixelement squared of Eq. ^ can be obtained from the expressions for the direct 
box of Eqs. (IB5I) . (IB6P by performing the following replacements: 

z ^ — z 

Mi _ Mi 

n.2 —<■ n.2 

Af' ^ Af' (B9) 
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and the four point functions Di and Dim are evaluated at [Di,D, 



Im 



[A, Dim]{mt, 0, 0, ml u, s, m|, m|, m\,mi). 

Appendix C: ANALYTIC EXPRESSIONS FOR THE NLO SQCD CORRECTIONS 

TO gg ^ tt 



The Feynman-diagrams of Figs. [5] and [6] represent the SQCD 0{as) contributions to the 
^luon fusion subprocess consisting of the following contributions: 

• the modification of the gti-vevtex in the s- and t(M)-production channels described 
by the UV finite (after renormalization) form factors Fv,Fm,Ga and p]''^^'^\ a^'^^''^\ 
respectively, 



• the modification of the ggg vertex {p^'^) and the gluon self-energy (11) in the s pro- 
duction channel, 

• the self-energy insertion to the off-shell top quark propagator in the t(u)-production 
channel (pf erf '^*'"'*), 

• and the box diagrams of the t and u production channels, (pf(* ^ c)' /'S^^**' ^7(1 b c)) 
(P°llb)^ ^!!(a,fe))' respectively. 

The form factors for the gti vertex, Fy, Fm, Ga, and the subtracted gluon self-energy, 11, are 
provided in Appendix B and Section III Bl respectively. The contribution of the renormalized 
ggg vertex of Fig. [29] to the NLO SQCD matrix element of Eq. ([5]) reads as follows: 



1 



q — u,c,t j 
d,s,b 



\ 4:Ul~ 

— (§ + 8m? )Eo(s, m? , m? ) ^5o(0, m? , m? ^ 



+ ^% Co (0, s, 0, m?^, , m| , m\ ) + 



3s 



(4m? - s)5o(s,m?,m?) 



4m 



^Eo(0,m?,m?) + m?Co(0,s,0,m?,m?,m?) + 



5Zi 



-gi -g, ■ - g - ^ \ - T 1 - J ■ ■ - gi -gj -g/ ■ -j^^ 

with 5Zi of Eq. ( fT3l l and the Bq and Co functions in the convention of Ref. 0|. 



(CI) 
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9 



4^ 



9 



4^ 




Figure 29: Renormalized gluon-gluon-gluon vertex correction at NLO SQCD. Graphs containing 
squarks are summed over the squark-mass eigenstates j=L,R(l,2) (no(with) mixing) and the quark 
flavors q={u,d,s,c,h,t}. 



v,t 



1. Parity conserving coefficients to the SMEs ' 



The t-channel parity conserving coefficients to the SMEs M/^'* read as follows: 
Vertex corrections: 

The vertex corrections consist of the stop-stop-gluino, p^f^g, the gluino-gluino-stop, Pi,ggt, 
contributions and the counterterm, pi^cT, as shown in Fig. [27l 



V,t _ V,t . V,t , V,t 
Pi — Pi,ttg ' Pi,ggt ' Pi,CT 



(C2) 



with 



v,t 

Pl,ttg 



V,t 
P4,ttg 



V,t 
Pll,ttg 

V,t 
Pu,ttg 



V,t 
Pl6,ttg 



i=i,2 



J2 iK^'^t-i) (C2 + C22 + C7i2) 



J=l,2 



V,t 
'Pl,ttg 



V I IxjiCn + C22 + 2Cu + Ci + C2) - Xj^iCo + Ci + C2)] 



-4p 



V,t 

14,ttg 



with [Ci,Cin,] = [Cz,C/m](m2,0,t,m? m? m? ) of Ref . [601, and 



(C3) 



v,t 
P^,ggt 



v,t 

Pll,ggt 



i=l,2 



2Xj m-gmtCo] 

5^|A+K-t) [C2 + C22 + C: 



i=i,2 



- H + ^DC'o + it - ml)Ci + 2Coo - 5o(0, m|, m|) 



i=i,2 
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rrit 



.2) + At!^(Ci + C2; 



= -4pr4t (C4) 

with [Ci,Cim] = [C*/, C/m]("^^, 0, t, m? , m?, m?) and the counterterm 

^Pi.CT = -^PilcT = ^^Zv (C5) 

with 6Zv of Eq. (fTSl) . The parameters are given in Eq. (IA2p . 
Top quark self-energy insertion: 

The self-energy insertion in the off-shell top-quark propagator as shown in Fig. [5] with the 
SQCD one-loop corrections shown in Fig. [8] is described by: 

^Pi'' = -(t + m2)(S^(t) + 6Zv) - 2m2(S5(t) - SZy - S^K) - Sv^(m,^)) 
^p^i' = 2£(Sv(t) + 5Zv) 

+ ii + mf)iJ:sii)-5Zv-^s{m^t)-^vim^t)) , (C6) 

where the SQCD one-loop contribution to the vector and scalar parts of the top quark self- 
energy, Sy 5, and 6Zv are provided in Section Hi B[ 
Box contribution: 

The box diagrams of Fig. [6] can also be parametrized in terms of coefficients to the SMEs 
Mj^'* as follows: 
Box a: g - g ~ g - ij 

Pii = 1^0 + (m| - m't - ml) Do + {ml - m',) - 6D00 

- i (2^22 + D222) - 2 (t + ml) (D12 + D122) - 6D002 - 2mlDii2 + {s - 2m]) D123] 
+ At [-2mgmtDQ\] 

P2°f = - E 2^.^ [^00 + ^002] 
i=i,2 

P^i = + 5Z [("^t - (Di + /^2) + 2Doo + (2Z}n + 2D,3 + /^iii) 
+ (t + m^) D22 + (i + 5m]) D^^ + 2(I^ooi + 2^002) 

+ £^222 + {i + 3m]) Du2 + (3t + 2m,2) D122 + (Sm] - s) Du3 + (i - s + 3m]) D123] 
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= +Y1 '^K [^22 + 2Di2 + D222 + 2D112 + 4Di22 + 2D123] 
j=l,2 

Pn'l = - E {^^'^ [^0 + (m| - m|) Do + (t + 2m? - m^) D, 
+ (t + ml) {Du + /^i3 + ^112 + ^123) 



t {D2 + 2Di2 + 1^122) - 2^00 + QDooi + mf Dm - (s - ^mf) D^] 
nit 



- XJ \{m\ + ml- ml) Do - 2(t + mf) Di - 21 D2 - Co 



P?2a = + 5Z 4 |a+ [Z^oO + 2^001 + D002] + AT ^ D, 

3=1,2 ^ 

P?i = + 5Z 2 {A+ [2Di + D2 + 2Dn + D22 + 4/^12 + 2D13] 
i=i,2 



777 - 

+ A- ^[2^1 + ^2] 
^ mi 



PFm = - 5Z 4 {A+ [4Di + 2/^2 + 3^22 + 12^12 + 2Dm + ^^222 

+ 6 (Z^ii + 1^13 + /^112 + /^113 + ^122 + Di2-i)] 

Tfl ~ 1 

+ A" ^ [4Di + 2/^2 + 2Dn + D22 + 4Di2 + 2Dr^ \ (C7) 
^ mi J 

where the three- and four-point functions are denoted by Co = Co(0, 0, s, m?, m?, m?) and 

[A, A»n, Amn] = [Di, Dim, Dimn]{mf , 0, , mf ,i , s, m^ , mj, m"!, m^) , respectively. 

Box b: ij — ij — ij — g 



i=i,2 



i=i,2 



= + E 2\^[^2 + ^222 + 2 (2Di2 + D22 + 1^112 + ^123 + 2/^122)] 



i=i,2 



P?26 = - E 4 |-A;[2Aoi + D002] + A-^Doo| 

,=1,2 J 

pS*6 = - E ^ {^^[^2 + A22 + 2(Di + 2Du + 4Di2 + 2Di3 + D22 + i^iii 
i=i,2 

+ 3^112 + 3Di23 + 3/^113 + 3D122)] 

- \]^[Do + D22 + 2(2Di + D2 + Dn + 2Di2 + D^)] \ , (C8) 
mi J 
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with the 4-point functions [A,Am,Amn] = [Di,Di„„Dimn]{mf, 0, 0, mf, i, s, m?, m|, m|, m|). 
Box c: ij — ij — g — g and g — g — ij — ij 



i=i,2 



P?:' = [2^00 + 2^002 + 2^003] 



i=i,2 

- £(Di2 + 1^22) - 2D001 - 2/^002 - m'^iDni + 2Dn2 - 2Du3 - Dns + Du2 - D 



223 



- s{Di23 + Dn3) - t{Dn2 + 2D123 + A13 + ^222 + 2D122 + D223) 
P°c = + Y1 ^^J [^12 + D22 + Du2 + 2Du3 + + D222 + 2D122 + D223] 

P?i,c = + E 4\^^00 

i=i,2 

P?2,c = + E 4 [2/^00 + 2/^001 + 2/^002] + At!^2Doo| 

j=l,2 J 

= + E 2 [^1 + ^2 + Al + 2Di2 + D22] + A7^(/^l + /^2)| 

i=l,2 ^ 

pS?c = - E S + ^2 + 2Du + 4^12 + 2^22 + Dui + 3^112 + ^222 + 3/^122] 



i=i,2 



?77 ~ 1 

+ A-^ [Di + D2 + Ai + 2/^12 + 1^22] > . 
nit J 

with [A,Am,Amn] = [Di,Dim,Dimn\im^ , 0, m2, 0, t, u, m? , m? m? m? ). 



(C9) 



Triangle-box: L —tj — g 



u^,u ^ ^ /A;(Ci + C2)-AT!^Co|(m?,5,m?,m|,m^„mf.) . (CIO) 
,_i 9 L nit ) 



3=1,2 



2. Parity violating coefficients to the SMEs Mf'^ 

The t-channel parity violating coefficients to the SMEs M/'* read as follows: 
Vertex corrections: 



V,t V,t I V,t I V,t 

with 

i=i,2 
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a 



v,t 

U,ttg 



-V,t 
14,«g 



i=i,2 



(C12) 



with [Ci,Ci„^] = [Cz,C/m](m^^O,t,m?,m|,m|) and 



i=i,2 



-J ''J ' 



cr 



y,t 



i=i,2 



0" 



U,ggt 



3 

9-^/ [^11 ~ ^22 + Ci — 



i=i,2 



with [Ci,Cim] = [C/, C/m]("^^, 0, t, m? , m?, m?) and the counterterm 



2SZ, 



(C13) 



(C14) 



where the renormalization constant, SZ^, is given in Eq. ( fTSi l. The parameters Aj^ are given 
in Eq. (1A2]). 

Top quark self-energy insertion: 



(t-m,2)(S^(£) + 5Z, 



with the axial vector part of the top quark self-energy, T,a, of Eq. ffT4 
Box contribution: 
Box a: g-g-g-ij 

^i°a = [Co + {ml+ml-ml)Do + {ml + ml)D2-QDo^ 

J=l,2 

- t {2D22 + /^222) + 4m,2 Di - 2{i + m^) (D^ + D122) - 6 /^oo2 - 2 m^D 
+ {s-2ml)Di2z] 
-2 J2 [Doo + D 



(C15) 



112 



2,a 



^002 



i=l,2 



cr 



4, a 



5^ 2Af [2/^00 + (2(Dn + /^is) + Dm] 



J=l,2 
,2 



cr 



a 



6, a 
14,a 



+ (m^ - m?) (Di + D2) + {i + 5m2) D12 + 2Dooi + 4L'oo2 + ^ ^222 + (t + Sm^) D 
+ (3t + 2mf) D122 + (3m2 - s) Du3 + (t - s + 3m2) D123 + {i + rn?) ^22] 
= [^22 + 2/^12 + /^222 + 2Du2 + 4/^122 + 2/^123] 



3=1,2 



J2 2A/ [/^2 + 2Du + D22] 



(C16) 



j=l,2 
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with the three- and four-point functions are denoted by Co = Co(0, 0, s, m?, m?, m?) and 
[A, Am, Amn] = [Di, Dirn, Dirnn]{m'^ , 0, , J , s, m"!^ , mj, mj, mj) , respectively. 
Box b: ij — ij — ij — g 

^4°^ = Yl 4A,^[/^oo + 2^001 + ^002] 

^6°^ = Yl + A22 + 2{2Du + A2 + Dn2 + D^s + 20,22)] (C17) 

with the 4-point functions [A,Am,Amn] = [Di,Di„„Dimn]{mf, 0, 0, , i, s, mj, m|, m|, m|). 
Box c: g — g — ij — ij and ij — ij — g — g 

(^2°;* = - XI 4A/[I^oo + ^002 + ^003] 

j=l,2 

= X 2Af [2Ao + 2Aoi + 2Ao2 - ml{D, + D^) 

+ m?(Di + A + 2/^11 + 3/^12 + A2 + Dui + 2Du2 - 2Du3 - Dus + D122 - A23) 

+ i{Du + A2 + 1^112 + 2Du3 + Du3 + A22 + 2^122 + A23) + s{Di23 + ^113)] 

^6°c = Y 4A/[^i2 + A2 + Dn2 + 2Du3 + Dn3 + A22 + 2D122 + ^^223] 
i=i,2 

= J2 2A/[^i + D2 + Dn + 2Di2 + 2Di3 + A2 + 2D23] (C18) 

i=i,2 

with [A,Am,Amn] = [Di,Dim,Diran]{m^ , 0, , 0, t, u, m? , mj, m?, m? ). 

Finally, the w-channel contributions p[^'^^'^, pf^l'^^ and crf^'^'''", o"-^^^'^^ can be obtained 
from the t-channel form factors by replacing i with u. 

3. The color factors 

The summation (average) over the color degrees of freedom in the course of the derivation 
of the parton cross sections leads to the following color factors: 

c^(l) = E(/.,.i;^^) UaMT^iT = 3^5,, Tri^T'^} = 12 

a.i'.c; c,d 

3,1 

a,b,c; 
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a,b,c; 




c*(2) = 


3,1 

" ^ ] {^'^jmJ'ml) i~'^'^jm'^ml) ~ ^ (3) = 

a,b,c; 


lb 

" y 


c«(2) = 


~ ^ ] {^'^jmJ'ml) i~'^'^jm'^ml) ~ ^ (3) = 
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" ~3 
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3 1 
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:)* = pil) 


4(2) = 


Q 1 

-c*(2)-- J](.<5,,5,,) {-iTf^Tl,r = 
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Q 

r'P) + 1 


4(3) = 


3,1 

^^*(3)-^$^(^5a.5.,) (-^T,t„r;j* = 
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3,1 


|c'(3) + 1 


4(1) = 

4(2) = 


ic'd) 

ic'(2) + i 5;(tf„j^y) (-ffj^r„'„)- = 


>)-! 


4(3) = 


3,1 

^^'(3) + i$^(z5a.5.,) (-zT,t„T^,r = 

a,o,c; 


i.'(3) - 1 


4(1) = 

4(2) = 


i,* 



a,6,c; 


■ 1 


4(1) = 

4(2) = 


-c"(l) 
2 ^ ' 

Q 1 > 

-c«(2)--5^(z5,,5,,) (-z7;':„Tlr = 


|c"(2) + 1 


4(3) = 


^^"(3)-^$^(z5..5.,) {-^T^M* = 

a,b,c; 


^c"(3) + 1 


4(1) = 

4(2) = 


a,b,c; 
3,1 


^c-(2) - 1 


4(3) = 


^^"(3) + ^5^(^5.,5,,) (-^Tf^T^,r = 


^c"(3) - 1 



(C19) 



(C20) 



a,b,c; 
3,1 
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